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We present the first multi-fluid analysis of a dense neutron star core with a deconfined colour- 
flavour-locked superconducting quark component. Accounting only for the condensate and (finite 
temperature) phonons, we make progress by taking over results for superfiuid ''He. The resultant 
^--s ■ two-fluid model accounts for a number of additional viscosity coefficients (compared to the Navier- 

^\] ' Stokes equations) and we show how they enter the dissipation analysis for an oscillating star. We 

, , provide simple estimates for the gravitational-wave driven r-mode instability, demonstrating that 

J^. the various phonon processes that we consider are not effective damping agents. Even though the 

results are likely of little direct astrophysical importance (since we consider an overly simplistic 
stellar model) our analysis represents significant technical progress, laying the foundation for more 
detailed numerical studies and preparing the ground for the inclusion of additional aspects (in 
particular associated with kaons) of the problem. 
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I. INTRODUCTION 



r^ 

f-^ I The state of matter at extreme densities continues to be an issue of vigorous investigation. The problem is conipli- 

I ■ cated, not only from the theoretical point of view, but also by the fact that laboratory experiments are restricted. For 

O ' example, while colliders like RHIC at Brookhaven, GSI in Darmstadt and the LHC at CERN probe hot quark-gluon 

*-J , plasmas they will never be able to explore the cold, extreme high density, region of the QCD phase diagram. In order 

^ ■ to test our understanding of the relevant physics we need to turn to astrophysics, and the dynamics of compact stars. 

—I, In fact, "neutron stars" represent unique laboratories of such extreme physics. With core densities reaching about 

one order of magnitude beyond nuclear saturation, they arc likely to contain exotic states of matter like hypcron 

phases with net strangeness and/or deconfined quarks. It is well-established that these states of matter should exhibit 

superfluidity/superconductivity at the relevant temperatures (neutron stars are born with temperatures ~ 10^^ K 

Cq ' and rapidly cool below ^ 10^ K). Moreover, observed radio pulsar glitches provide strong evidence for the presence 

of a partially decoupled superfiuid component in these systems. The modelling of the dynamics of these, potentially 

very complex, objects presents a serious challenge. 



A key aspect of the problem concerns the fact that a superfiuid system has additional dynamical degrees of freedom. 
f*~^ ' This is well-known from experiments on laboratory systems like ^He, which exhibit a second sound associated with 
(^ [ thermal waves [l|, |2| . Analogous "superfiuid" modes have been studied in detail for superfiuid ncutron-proton-clcctron 
mixtures relevant for the outer core of a neutron star [3-[3 ■ For simplicity, these studies have almost exclusively ignored 
thermal effects (the work in @ is a notable exception). While this is a useful first approximation it is clear that the 
zero temperature assumption must be relaxed in a realistic model. Basically, due to the density dependence of the 
various superfiuid pairing gaps (see [7[ for a guide to the literature), there will always be regions in a neutron star 
where thermal effects are important (in the vicinity of the critical density at which the phase transition occurs). 
Understanding the nature of these transition regions, and their effect on various aspects of neutron star dynamics, is 
one of the main challenges for research in this area. 

In this paper, we describe a first attempt at modelling thermal dynamics in a superfiuid neutron star. We focus on 
stars with a colour-flavour-locked (CFL) superconducting quark core Q at finite temperatures. This is an interesting 
problem for several reasons. First of all, the simplest possible model for this system considers a quark condensate 
coupled to a gas of phonons. This problem is analogous to ^He at low temperatures, and hence we can bring our 
recent dissipative two- fluid model Q to bear on it (more or less directly). The lessons we learn from this exercise 
should inform the development of finite temperature models for the superfluids in the outer core and the neutron star 
crust. Secondly, even though there have been discussions of the dynamics of the different CFL phases, in particular 
in the context of the gravitational- wave driven r-mode instability (see [lO| for references), the supcrfluid aspects have 
(as far as we are aware) not previously been accounted for. Our discussion begins to address the relevance of the 
additional degrees of freedom in these systems, and provides some insight into the nature of the different "fluids" 
involved. 

As an application with immediate astrophysical relevance, we will work out the inertial r-modes and estimate the 
relevant viscous damping rates for the "simplest" model of CFL matter. We consider the "cool" regime where the 
temperature is significantly below all the quasiparticle energy gaps. In this regime, dissipation may mainly occur 
due to phonon interactions. One reason for considering this model is that there are results in the literature for 
both bulk- and shear viscosity [ll|, [l3 as well as the mutual friction associated with supcrfluid vortices [l^- Of 



course, the simple "condensate plus phonon" model that we consider is not the whole story. It should apply at 
asymptotically high densities, but may not be the true ground state at lower densities. The discussion in p^| - [l7| 
adds extra dimensions to the problem by considering the bulk viscosity due to kaons (allowing for flavour-changing 
processes), which will be present at higher temperatures. At first sight, this mechanism will only be relevant for very 
hot stars (T > 1 MeV ~ 10^° K) since it assumes that there is a thermal population of kaons. Below the critical 
energy where the kaons appear the contribution to the bulk viscosity is exponentially suppressed and may not be 
that important. However, in more recent work on the so-called CFL-K*^ phase [16l |. it is argued that the main low 
temperature mechanism involves condensed kaons. This is important since the kaon condensate will remain present 
as T — > 0. The upshot of this is that the problem requires a "multi-fluid" analysis at all temperatures. The simplest 
model would have three components; the quark condensate, the kaon condensate and finite temperature excitations 
(phonons and thermal kaons). This problem is more involved than the case that we focus on here. Nevertheless, our 
results provide an essential starting point for investigations into the dynamical role of the kaons. Most importantly, 
by providing the "hydrodynamics" view of the problem we illustrate the input needed to study various dynamical 
scenarios. This should stimulate further discussion between experts on different aspects of this multi-faceted problem, 
as required to make progress in the future. 

II. FLUX-CONSERVATIVE TWO-FLUID MODEL 

We consider a deconfined quark system system that contains a single CFL superconducting condensate and a 
phonon gas. Formally, this problem is identical to that for ^He at low temperatures. Hence, we can take as our 
starting point the recent discussion of superfiuid Helium [9| , which is based on the convective variational approach 
to multi-fluid dynamics [18|-|20[ and which we know is in one-to-one correspondence with the orthodox formulation 
developed by, in particular, Khalatnikov [ij. 

The variational model takes as its starting point an energy functional E, representing the equation of state. This 
energy determines the relation between the various fluxes and the associated canonical momenta. In the present 
case, we distinguish between the massive particles in the system, with number density n and flowing with nvf from 
a massless entropy component with number density s and flux sv^. The latter represents the phonons (which are 
treated in the fluid approximation, i.e., we assume that there exists a suitable defined "avera ge" transport velocity 
for each of the two components in the system) The two momentum densities then follow from 51 1 
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and 



< = 2awr , 



(2) 



where wf^ = vf — wf and a represents the entropy entrainment 2l|. It has been assumed that E = E{n, s, w'^), as 
expected in an isotropic system, and we have defined 



dE 



dvj^ 



(3) 



omitting the indices on w for clarity. 

As discussed in [9|, Il9| , the associated momentum equations can be written 



/f = at< + V,(« + D%) + nV. (^Mn - \mv^}j + 7r;^V,< , 



and 



(4) 



where fi^ is the matter (quark) chemical potential; 

'dE 



Hn 



dn 



3 

s ' 



(5) 



(6) 



and we have used the fact that the temperature corresponds to the entropy chemical potential |2l| 



(7) 



In these expressions, Df, represent the viscous stresses while the "forces" /f ahow for momentum transfer between 
the two components. In the foUowing we wiU assume that the system is isolated, which means that /" + ff = 0. As 
we will discuss later, the force terms can also be used to account for "external" forces like gravity. 
In the present context, when we are dealing with a single matter quantity, we will have [19| 

dtn + \/j{nvi) = T^^O, (8) 

as there is no particle creation/destruction. At the same time the entropy can increase, so we have 

dts + yj{svi) = r,>o . (9) 

We also have [l^ 

TV, = -f>L - D{,V,vl - mv,ui. , (10) 

where 

D,,^D^^ + iyi^. (11) 

For an isolated system we know that, if we impose the superfluid constraint of irrotationality the number of 
dissipation coefScients reduces significantly. As discussed in Q , we have 

- A, = 9^, {C^if + Ce.) + 2i^Q% . (12) 

where we have defined j* — nwl^^. We also have; 

i (/f - \/iD^\) = v.* (13) 

with 

^ = C""Vu-' + C"0s . (14) 

At this point, only four dissipation coefficients remain in the problem. 
In the above expressions, we have used the standard decomposition; 



S/,v';=&l^ + -g.,Q. + e,,kW,' (15) 



in terms of, the expansion 



the trace-free shear 



and the "vorticity" 



e, ^ v,vi , (16) 



% = U^^^'j+^M-l9uQs] , (17) 



W: = ^e^^H^jvl^VkV^^). (18) 

We use analogous expressions for gradients of the relative velocity. The definition of the various quantities should be 
obvious from the constituent indices. 

However, since the system that we consider is not irrotational we need to consider relaxing the assumptions on the 
dissipation coefficients. This involves making some subtle decisions. In the case of an irrotational flow it is natural to 
assume that /„ = 0. When the superfluid rotates, and vortices are present, the force will not vanish. It is necessary to 
account for dissipation due to, for example, the scattering of phonons off of the vortex cores. The standard approach to 
the rotating problem is to add in this "mutual friction" force, keeping the other dissipative terms as in the irrotational 
case. This strategy ignores a number of dissipative terms that would, at least in principle, be allowed in the equations 
of motion 13,[i^. The role of these additional terms has not yet, as far as we are aware, been investigated. 



Ignoring the potential relevance of most of the additional dissipation channels in the irrotational case, we simply 
account for the presence of vortices by i) assuming that ^ represents fluid elements with both a condensate and a 
smooth-averaged vorticity arising from the vortices, and ii) accounting for the vortex mediated mutual friction by 
allowing a force [22| 

/r' - B'pnn.ue,,kii^wt + Bpnn,€,,ke'''"'k^Kiw^t , (19) 

to act on the particles (with a balancing force affecting the excitations). Here n„ is the vortex area density and the 
vector K* = KK* (the hat represents a unit vector) is aligned with the rotation axis and has magnitude k = h/2m. 
Since we will only consider the effect of phonon scattering off of vortices, we expect to be in the weak mutual friction 
regime where B' <C B. This means that the first term in the force ([T^ can be ignored, leaving only the second, 
dissipative, contribution. 

III. ROTATING EQUILIBRIUM MODELS 

In order to set the stage for the discussion of r-modes, we need to provide a suitable rotating equilibrium configura- 
tion. To do this, we note that the two components flow together (there is not heat flux) when the system is in thermal 
equilibrium. This means that wf^ — 0, which implies that irf = 0. Thus, it follows from the entropy momentum 
equation ([5]) that wc must have 

sV,T = — > T = constant . (20) 

This is trivial; if the system is isothermal then there will be no heat flux. 

Rewriting the remaining momentum equation (j3]) , making use of the continuity equation, we find (after accounting 
for the gravitational force as the gradient of the gravitational potential $) 

dtvf + viy.vf + V,(/i„ + $) = , (21) 

where /in = P-n/^n- Using the Lie-derivative along the fiow, £„^ , we can rewrite this equation as 

dtvl + £,„< + V, h"i„ + $ - ii)2 j = . (22) 

Restricting ourselves to stationary models we have 

dtv^ = . (23) 

Next, uniform rotation implies 

< = f^ne^ , (24) 

and, for axisymmetric models, it then follows that 

C^vf = . (25) 

Hence, the required equilibrium models are determined from the usual Bernoulli-type equation 

/in + $ - -Wn = constant . (26) 

This analysis shows that a rotating configuration can be obtained in the usual way. Once we assume thermal 
equilibrium, we are dealing with a single-fluid problem. 

IV. LINEAR PERTURBATIONS 

Our main aim is to develop the tools required to make quantitative estimates for the r-mode instability in neutron 
stars with a CFL core. In doing this, we want to account for the fact that CFL matter requires a multi-fluid 
description. The analysis proceeds in three steps. First we need to formulate the linear perturbation problem for the 
system. This step is interesting because, as far as we are aware, this is the first time that perturbations of a matter 



plus masslcss entropy system have been considered in an astrophysical context. The results provide useful insights into 
finite temperature superfluid dynamics, and should be relevant in a broader context. The second step corresponds to 
determining the pulsation modes of the system, in this case the r-modes, and finally we need to estimate the damping 
timescales associated with the different dissipation channels. In the discussion below, we only work out the first of 
these steps in detail. In order to obtain useful estimates for the r-mode instability we simplify last two steps by 
considering a uniform density model. This has the advantage that the r-mode solution is simple, and we can evaluate 
the dissipation rates analytically. In a sense, we do not expect this approximation to be too bad because it is well 
known that the density profile for a canonical 1AMq/10 km strange star (described by the MIT bag model) is almost 
flat (see [H). 

Of course, we really need to develop a consistent model for a realistic equation of state. This is essential if we want 
to study hybrid stars, where the CFL phase is only present in the core. However, the problem is complicated by the 
fact that the bulk viscosity damping of the r-modes requires the perturbations to be worked out to second order in 
the slow-rotation approximation |24l |. This necessitates a numerical solution, which makes some of the qualitative 
behaviour of the results less clear. Another complicating factor is the well-known fact that it only makes sense 
to use a realistic equation of state in a fully general relativistic analysis (see |25| for discussion). The r-modes in 
single component relativistic stars have been studied |26l - l29| , but there has not yet been any serious analysis of the 
corresponding multi-fluid problem (although see |30|). A focussed effort in this direction should be encouraged. 

A. The non-dissipative problem 

In order to work out the r-mode solutions, it is advantageous to work in a rotating frame. Focussing, for the 
moment, on the non-dissipative equations we find that, in a frame rotating uniformly with r2* we have 

dtirf + Vj («) + nW, f Mn + m$ - ^-mvA + 7r"V,< + 2pe,,kn^v^; = , (27) 
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dtn! + Vj (viTrl) + sV,T + Tr^^V^v^ = , (28) 



and 



It is notable that the Coriolis force does not affect the entropy equation (|28|) . This is, however, not surprising. It is 
well-known that incrtial forces are proportional to the mass, and since our entropy component is taken to be massless 
it should not be affected by the Coriolis force. From a technical point of view, it means that the problem we consider 
is subtly different from the two-fluid r-mode problem discussed in [3l|, (sj . 

Let us now consider perturbations of the rotating equilibrium models discussed in the previous section. Considering 
Eulerian perturbations (denoted by S) we have, in the rotating frame (using the Cowling approximation (5$ = 0, and 
not explicity denoting the velocities as perturbations since they vanish in the background anyway); 

pdtv^ - 2adtwT + py^6^1n + 2pe,jk^^vl = , (29) 

and 

2adtwT + sV,5T = . (30) 

Adding these we get an equation for the total perturbed momentum 

p9t< + pV,5jl,, + sV,5T + 2pe,,k^^vl = . (31) 

Noting that the pressure p is defined by (for a co-rotating equilibrium model) 

V,p = pV./in + sV,T , (32) 

we see that 

Vi5p = 5pVi/i„ 4- pViSpr, + SsViT + sV,6T = — V,,p + pV^Spn + sV,5T , (33) 

P 

since V^T = in the background. Hence, we have the usual Euler equation 

pdtv'l + V,5p - ^V,p + 2pe,jk^^vl = . (34) 



We also have the perturbed continuity equation for the particles 

9t(5p + V,(p<) = 0, (35) 

and a conservation law for the entropy 

dtSs + V,(s«:) = dtSs + V,(s<) - V,(s<J = . (36) 

The entropy is conserved since there is no heat flux in the background (Fg is quadratic in the heat flux [2l|). 

Once we provide an equation of state for matter, we have all relations we need to study the linear dynamics of a 
non-dissipative finite temperature CFL quark core. 

B. Energy integrals and dissipation 

Our main aim is to work out the r-modcs and establish the parameter range in which gravitational-wave emission 
triggers a secular instability. In order to assess the relevance of this instability we need to consider the various 
dissipative mechanisms that counteract the growth of an unstable mode. The damping due to shear- and bulk 
viscosity in a two-component system can be estimated using the strategy set out in J3ll |32|. That is, we use energy 
integrals to estimate the various timescales. In essence, the damping timescale r associated with any given process is 
obtained from 



2E 



dE 



dt 



(37) 



where E is the energy associated with the flow and dE / dt is the rate of energy loss due to dissipation. This estimate 
should be accurate as long as the dissipation is weak enough that it docs not affect the flow on a dynamical timescale. 
The energy associated with a given perturbation can be obtained from the equations of motion (|29[) and (|30l) . 
Multiplying the first equation by ujj and the second by vl (where the bars denote complex conjugates) and adding, 
we find that the kinetic energy is given by 



i^=i /p(kP-yhnsnd^. (38) 



2 



In the case of the r-modes, this is the leading order contribution to the energy. In fact, it turns out that the second 
term in the bracket is of higher order in the slow-rotation approximation (for the same reasons as in j32l|). Hence, the 
r-mode energy is well approximated by 



E^\J p\v^\^dV . (39) 



In order to evaluate the energy loss due to viscous damping we add the dissipative part of the stress tensor Dij to 
the equations of motion. The total energy dissipation then follows from the body integral of 



-viy^Dt,~^'iyi, 



After some algebra, this leads to 

dE__ 

'dt ~ 

where we have defined 



I [vQi^eij + C"" |e|' + 2CRc (eOs) + c\&s\^]dv , (40) 



e = v,/ . (41) 



Recall the relations (fT6)) and (|17p and the definition j' = nwl^^. 

In order to combine this result with existing results for the viscosity coefficients we need to translate our variables 
into those of the "orthodox" two-fluid model due to, for example, Khalatnikov [l|. The required translation has been 
discussed at length in [9| , and the key relations that we need are given in the Appendix. Carrying out the comparison, 
we find that the shear viscosity terms are exactly the same in the two descriptions. Hence, we can use rj from [IJI 



without change. In the case of the bulk viscosity, wc find that the two sets of coefficients have different mass scahngs, 
and we identify 

C = C2 , C = mCi = mC4 , C"" = m^G • (42) 

This means that, when expressed in terms of the standard coefficients, the bulk viscosity dissipation rate becomes 

dE 



dt 



J [rjei^ei^ + m^Cs |e|' + 2mCiRe (96,) + C2 |e,|'] dV . (43) 



To conclude, there are three bulk viscosity coefficients rather than the usual single one, and both dynamical degrees 
of freedom are needed if we want to evaluate the dissipation integrals. Note that the dissipation is expressed in terms 
of J* and vl, rather than wj, and u;,\g, the variables that we will solve for when we determine the r- modes. It is, of 
course, straightforward to express (|33]) in terms of these variables. Unfortunately, the dissipation integrand is then 
rather messy. We find 

+ C2 |e„|' + [C2 - 2pCi + p'Ca] |e„/ - 2(C2 - pCi)Re (enOn,) 

+ 2CiRe [On «V,p)] - 2(Ci - PCs) [e„. «,Vjp)] + Cs Kv,p|' }dV . (44) 

where 

e„ = Vj< and e„, = Vjwi, . (45) 

However, as we will see later, this expression simplifies considerably in the r-mode problem. 

Before moving on, it is worth making a general point. It might be tempting to suggest that the terms involving 
wjjg in (j44p should be less relevant than those involving only v"^. However, without solving for an actual oscillation 
mode one cannot make this argument precise. In general, the bulk viscosity contributions from terms involving wl^^ 
cannot be neglected. Since a mode-oscillation typically involves both degrees of freedom (unless the star is unstratified 
[33I, |3Jj), we need to determine the nature of the oscillations before we make further simplifications. 

V. A SIMPLE MODEL EQUATION OF STATE 

In order to obtain quantitative results for the r-mode instability we need to provide an equation of state. As our 
main focus is on the two-fluid aspects of the problem (ignored in previous studies) , the equation of state must account 
for the thermal excitations, i.e. the phonons. Moreover, if we want to investigate the relative importance of the 
different bulk viscosity terms we need to allow for a relative flow at the perturbativc level. This is also important 
if we want to consider the superfluid mutual friction [13|. However, if we insist on the model being truly "realistic" 
then the problem becomes much more challenging. Hence, we will make a number of simplifying assumptions. This 
is not only practical, it is also quite reasonable since this is our first attempt at a quantitative analysis. Future work 
should aim to relax some of our assumptions. 

As discussed in Section II, the equation of state that we require takes the form of an energy functional E = 
E(n, s, w'^). Once this functional is provided, we can work out all the quantities that we need to construct a rotating 
model and calculate the r-modes. Unfortunately, we do not have an equation of state of the required form. In fact, 
some of the parameters that we need are usually not considered, essentially since they are not required if one is only 
interested in equilibrium configurations. In order to make progress we have to build a suitably simple model equation 
of state. To do this, we assume that the thermal contribution due to the phonons adds to a zero-temperature equation 
of state for the condensate. This strategy has previously been developed for a relativistic superfiuid [S^, and builds 
on the classic expressions for a non-relativistic phonon gas, see [l|. 

A. The phonon gas 

Let us assume that, at zero temperature matter is described by a barotropic model such that E{n, s — 0) — Eo{n). 
Then the chemical potential follows from 



the pressure is obtained from 

dpo ~ ndfx , (47) 

and it is easy to work out the speed of (first) sound; 

Co = ^ , (48) 

Now consider thermal excitations represented by a phonon gas with a hnear dispersion relation, with slope cp . 
Working out the thermodynamics of such a gas one can show that its contribution to the pressure, ip, is given by [35| 

^ 45 (2^fico)3 V cl) ■ ^^^^ 

Here, and in the following, it is to be understood that w'^ corresponds to if^g. From this we obtain the entropy density 
(in the matter frame) via 

This leads to the explicit result; 

_4^_167r5 / kT Y f w'Y^ 



ds 2tt^ fkT^^ 



which means that the heat capacity is given by 

in agreement with the phonon result in [l6|. 

Finally, one can show that |35| the "normal fluid" density, pN, required in the orthodox superfluid formalism (see 
Appendix) is given by 



45 {2nh)^ \^) V'^J ~o 



P^ = -^7KZ^J-] (^--a) -• (53) 



This result is consistent with 



4Ep„ /, w'-- ' 



^N^J^ 1-^1 . (54) 



\ "-0 

where the phonon energy is 

In our formulation of the problem, we need the entrainment between particles and entropy. As discussed in [9| the 
relevant entrainment parameter is related to pTss according to 

a = --[l--) «--, (56) 

where the last approximation is accurate in the low temperature limit, when p^ <^ p. This is the regime that we are 
considering here. 

It is worth noting that the above results lead to 



In order to be consistent we need pN "C p « 7 x lO^^g/cm (the average density for a canonical neutron star), which 
translates into T ^ 2 x 10^^ K. In other words, the model should be valid for all astrophysical neutron stars. 

Let us now use the phonon gas results to construct a "complete" model equation of state. That is, we want to 
deduce a consistent energy functional E. This functional should be such that the temperature is obtained from 



dE 
97 



Now, combining ([^^ with ([5T|) we find that 



where 



dE_ 



- 1,1/3 



B 



IGtt" 



5 \ 1/3 



1 



45 J 2nhco 
Since cq ~ cq^u) (by definition) we can integrate to get 



w 



-2/3 



E = Ei{n, w^) + -^s^/^ ^E,+3^ = Ei+ E^^ (l - ^ 



(58) 



(59) 



(60) 



(61) 



We learn that in the limit of a low relative velocity {w -^ cq) the equation of state is simply given by 

E = Ei+ Eph ■ (62) 

This result is quite intuitive. Using the fundamental relation of thermodynamics 

p + E^n^i + sT , (63) 

it is also straightforward to show that the total pressure is given by 

p = Pq + 11! , (64) 



as one might have expected. 

Now that we have the required energy functional, we can determine the entrainment parameter from 



dE 



dw"^ 



Working this out, we find that (after some algebra) 



Comparing to (j56p we see that the model is consistent with 

dEi 



dw"^ 







(65) 



(66) 



(67) 



Hence, it is natural to take Ei = Eo(n), i.e. simply add the thermal phonon energy to the zero-temperature equation 
of state. The final result is then 



(68) 



E^Eo + 3^^Eo + ^ph 1 - 



In order to proceed, we need to provide the zero-temperature equation of state. A "realistic" model should obviously 
be based on QCD. It is natural [l3| to use the MIT bag model. In the case of CFL matter, this leads to (using units 
where the speed of light is unity, c = 1) 



Po«g-i^ + g(4A^-™D 



(69) 
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where B is the bag constant and A is the pairing gap associated with the CFL condensate. 
For simplicity, we wiU use 

Po«i(p-4i?), (71) 



in which case the speed of sound is constant; 



t4- 



As discussed in [10l.ll7l|. this approximation may be quite good, essentially because the two contributions to the last 
term in (|69|) almost exactly cancel each other. Later, when we work out the various r-mode dissipation integrals, we 
will simplify the problem further by assuming that the density is uniform, with a constant speed of sound. 

B. Viscosity 

In order to determine the damping timescalc for the unstable r-modes, we need to supplement the non-dissipative 
model with viscous terms. From the discussion in section IIIB we know that we need one shear viscosity coefficient and 
three bulk viscosity coefficients. In order to account for dissipation associated with the superfluid vortices, we should 
also consider the mutual friction. Some of the required viscosity coefficients have been discussed in the literature. In 
particular, the dissipation due to various phonon interactions has been investigated jlll4l3l |. In fact, the availability 
of these results is one of the key reasons for us focussing on the condensate plus phonon model. Having said that, the 
available results are incomplete. To what extent this is the case, and how we have dealt with this problem, will be 
discussed below. 

Before we proceed it is worth iterating the point that we are focusing on the phonons because they represent the 
simplest in a hierarchy of relevant multi-fluid problems. The presence of kaons, either as thermal excitations or a 
condensate [1^ [l^ , would require us to extend the model to account for additional degrees of freedom. The issues 
involved are similar to those for a superfluid hyperon core |25| , and we expect to consider them in the future. 

Let us first consider the phonon shear viscosity. The relevant viscosity coefficient is determined in [l2|. We should 
use (following the discussion above, we are using Cq = 1/3 here) 

77 = 2.46x10^6/ ^h ^YH!^^ g/cms, (73) 

' V300McV/ V ^ / 

where /ig is the quark chemical potential. We arrive at this result by taking both /.Jn and the quark "mass" m in the 
fluid model to be equal to the quark chemical potential ^q. The temperature scaling of "q is the same as in the case of 
''He, which makes sense given that the involved phonon processes are the same. However, in the context of neutron 
stars, this model is severely limited. Basically, the "hydrodynamic" treatment of the phonons is no longer valid if 
their mean-free path is larger than (or comparable to) the size of the system. As in kinetic theory, the mean-free path 
associated with the shear viscosity follows from (up to a factor of order unity) 

?7 ~ pnCqA . (74) 

For the phonon gas model this leads to 

A.5xl0-(^L)Y1^V-. (75) 

V300MeV/ \ T J ^ ^ 

In other words, for a 10 km star the model will not be valid below ~ 10^" K. This is obviously a problem, since 
mature neutron stars are expected to be significantly colder than this. In essence, we ought to model the phonons as 
ballistic in the low temperature regime. However, there is some evidence from laboratory experiments on ^He that 
the fluid model remains relatively accurate also at lower temperatures. Given this, we will use the two-fluid model 
also in the regime where it is no longer formally valid. Still, in order to describe such systems we need a different 
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model for the shear viscosity. In order to make progress we wih simply assume that the mean free path is limited by 
the size of the system. That is, we will use an effective shear viscosity given by 

Vcs ^ PnCqR , (76) 

where R is the radius of the star. This may seem like a rather drastic assumption, but it has recently been shown [S^ 
that in the case of Helium it leads to results that agree quite well with the observed sound attenuation. In particular, 
we would have rjcs ^ pj^ ^ T'^ at low temperatures, which means that the viscosity weakens as the phonon density 
decreases. This makes (at least qualitative) sense. Explicitly, we get 

ryeff « 5 X IQi^ (^^^^ g/cms . (77) 

This model produces a maximum in the viscosity as a function of temperature, a feature that seems quite natural. 
Having said that, more detailed work on the low-temperature phonon problem is obviously needed in order to improve 
on our results. This is, in fact, a very interesting problem; when the mean-free path is large, the phonons interact 
with the "surface" more frequently than with each other. In principle, the hydrodynamics approach should not be 



valid in this regime. Yet, there is some evidence from studies of heat conduction in nano-systems (see 21| for a recent 
discussion) that a judicious choice of surface boundary condition for the phonons leads to a useful "fluid" model. We 
plan to explore this idea further in the future. 

We now turn to the bulk viscosity. Only one of the three required coefficients has, so far, been calculated in detail. 
In the static limit, the analysis in [11| leads to the result 

&(. = 0,.1.75xlO-(^5i^)'(i^)g/c™, (78) 

where rris is the strange quark mass. However, in order to use this result to study oscillation modes, we need consider 
the frequency dependence in more detail. This is essential since the bulk viscosity is a "resonant" mechanism that is 
particularly effective when the involved dynamics has a timescale similar to the relevant relaxation time. The static 
limit only provides partial information. Following |14{ we will use 

C»- ,^^, 2 «»(T), z = l-3, (79) 

where the "amplitudes" a^, in general, depend on the different physical scales of the system, while r is the relaxation 
time for the processes that give rise to bulk viscosity. According to [1J|, the relaxation time scales according to 



r - -7^ ■ (80) 



The scaling of all three viscosity coefficients, in the static limit, has been determined in [ij]. We should have 



Ci « a^-^ , (81) 

-t fJ-q 

C2 « 52^ , (82) 

C3 « 53-^ , (83) 

^ Pq 

where a^ are constants. This, together with the relations in ([75)1 and ([50)1 allows us to determine the frequency 
dependent bulk viscosity coefficients, up to the unknown amplitudes ai and /?. These can not be determined with the 
approximations used in |14l |. although we can obviously infer 6(2 from (|78p . From the analogous problem for ^He, we 
know that the three bulk viscosity coefficients may be of the same order of magnitude. In absence of detailed results 
for the quark case, it thus makes sense to assume that Ci and Cs have a similar amplitude to ^2- In particular, this 
would suggest that we parametrise Ca, in the static limit, according to 



10i4g/cmV \ T J \ M, 



C3(^ = 0) L.,/, 3 ) ="X10^'(^7^) ( ^) g/cms, (84) 



where a is an unspecified parameter. The frequency dependent result then takes the form 



C3 



1 + u;2t2 ' 



with 



= /3 X 10^ 



7 /lO^ K 



V300 MeV/ 



T / V 300 MeV 
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(85) 



(86) 



Here, (3 is another undetermined parameter. In principle, one would expect both a and /3 to be of order unity, but 
given the lack of precise information we can vary them and assess how this affects the bulk viscosity damping of 
the r-modes. Note that the coefficient in (j85p has a maximum at the resonance frequency w = 1/r. We know from 
previous work that it is essential to understand how this resonance frequency relates to the r-mode frequency. A precise 
statement to this effect is not possible, given the free parameters in our model, but the strong temperature dependence 
in (j86p means that even a change of several orders of magnitude in r would produce a relatively small change in the 
temperature at which the resonance occurs for a given mode frequency. This suggests that the uncertainty associated 
with /3 may not affect our analysis too severely (assuming that the temperature scaling is correct, of course). If we, for 
example, consider a mode frequency lu = 10^ s^^ and an amplitude such that r = 10® s, we would have a resonance 
around T w 3.5 x 10® K. If, on the other hand, we take a drastically shorter relaxation timescale of r = 1 s, then 
the resonance appears near T « 2.8 x 10^ K. That is, the resonance temperature would shift by less than an order 
of magnitude. This model is obviously phenomenological, and allows us to proceed, but an actual derivation of the 
frequency dependent viscosity coefficients is needed if we want more detailed results. 

It is obviously important to compare our results to other relevant estimates for the r-modes. The natural, and most 
immediate, comparison would be to unpaired quark matter. In that case, which corresponds to a single fluid problem, 
we have [11 [H 



Tj^lAx 10^ 



17/O.r 



5/3 . J, X -2 



where a^ is the strong interaction coupling constant, and 

A{T) 



j/cms , 



(87) 



C 



uj^ + B{T) ' 



where 



and 



A(r) = 1.8 X 10^^ ( ^^ Y( "^- Y(^] g/cms^ 
^ ' V300MeV/ VlOOMeV/ ^10 K/ ^' 



B{T) = 2.6 X 10^ 



Mq 



aOOMeV 



^K 



T 



109 K 



(89) 



(90) 



We also want to quantify the relevance of the vortex mutual friction for the r-mode instability in CFL matter. In 
order to do this we need a representation of the counter-moving degree of freedom associated with an r-mode (which 
is of order $7^ ^0. the slow-rotation approximation) [32|- Before we consider this problem, we need to translate the 
mutual friction results from [l3| to our formalism. In order to relate the parameter B to the results in [13, let us 
consider equation pop where we now include a mutual friction force of the form l\\9\ . This leads to an equation for 
the evolution of the relative velocity of form 



dwf 
~dt 



Za 2a 



JjfeK Wjjg 



+ Bn^e,jke'''^k^Kiw'^] 



(91) 



Given this equation, and the results in the App endix, we can compare our formalism to the standard mutual friction 
description, see e.g. equations (2.2)-(2.3) in [3^. This allows us to determine B from the parameter a that is calculated 
in 111. This leads to 



B 



2^ 
405 



2a 
P 



(1 



(92) 
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We see that the mutual friction vanishes as the temperature decreases. This obviously makes sense since there will be 
no phonon-vortex scattering when the phonon gas becomes dilute. Moreover, according to this result, we are safely 
in the extreme "weak drag" regime (c.f. the discussion in [331 ) where 

6' « 6^ < 1 . (94) 

This implies that the phonon mutual friction will not damp the r-modes efficiently (probably as expected) [32| . This 
result is confirmed by the detailed analysis in section FVIII However, these results come with an important caveat. 
Strictly speaking, our analysis is only valid as long as the phonons can be described as a fluid. As we have already 
explained, this is not the case for mature neutron stars. In fact, the results of [l3| are supposedly derived for ballistic 
phonons. This means that our analysis is somewhat inconsistent. At temperatures above (say) 10^" K the mutual 
friction parameter from jl3| may not apply, and below this temperature our two-fluid model may not be appropriate. 
However, the latter issue may not be that important. After all, one would expect the phonon mutual friction to 
weaken dramatically at lower temperature, c.f. (j93p . meaning the any quantitative errors in the analysis will be of no 
real importance. 

VI. APPROXIMATE R-MODE RESULTS, SHEAR- AND BULK VISCOSITY 

Our (simple) model for a two-component cool CFL core is now "complete" , and we may conduct the r-mode analysis 
as in [3l|, [S^] . In order to quantify the damping due to bulk viscosity, we need to account for terms of order fl^ in 
the analysis. This complicates the problem, since the centrifugal deformation of the star's shape enters at the same 
level. In view of this, and the fact that this is a flrst exploratory study, we will make a sequence of approximations 
that allow us to proceed analytically. This strategy also makes sense since we are using a simplified equation of state. 
The results of our analysis should help determine whether a full numerical study is worthwhile. 

First of all, we need to make a choice of primary variables. The linearised equations link the four scalar variables 
(5p, ST, Ss and Sp. In principle, we can use the equation of state to express any two of these in terms of the other 
two. Following the strategy set out in |32| we will work with 5p and ST. For our model equation of state we have 
(since the background is co-rotating) 

Sp = Spo + (^] ST = Spo + sST , (95) 



which means that we can express the equations in terms of Spo and ST. We then have 



^P^[^) ^Po + (^) 5T = \sp^ , (96) 



since cq is contant, and 



^-(^)*" + (§)^^^l^^' <^" 



since po =Po{n). 

In order to avoid numerics we will make use of a "trick" used by Lindblom, Owen and Morsink in one of the 
early r-mode instability papers J4C| . The basic idea is to neglect the rotational change in shape in different terms in 
the perturbation equations. Once this is done one can "estimate" the bulk viscosity from the leading order r-mode 
solution. Although this simplification is not not based on a rigorous argument, one can show that the bulk viscosity 
damping timescale is of the right order of magnitude (within about a factor of 5 of the true second-order slow-rotation 
result |4l|). This is good enough for our present purposes. 

In our case, the calculation would proceed as follows: First consider the continuity equation 

iujSp + pV,< + v\y,p = , (98) 
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where we have assumed that all perturbed quantities behave as e*('^*+™''') . If we simply omit the last term (that would 
vanish for a spherical star sinee an r-mode is purely toroidal to leading order), we get 

V,<«--<5p=-^<5po. (99) 

P pet 

This may not be a very accurate estimate, but it allows us to progress without having to determine the rotational 
corrections to the r-mode (which depend on the term that we have neglected). Now consider the second degree of 
freedom in a similar way. From the entropy equation we get 

iujSs + s (V,< - V,<,) + « - <jV,s = . (100) 

In fact, for the model equation of state the last term vanishes identically sinee s ~ s{T) and V^T = for our 
background model. It follows that 

5s 
V,< - V,<, « -iu:— . (101) 

s 

The approximation is completed by the i^-components of the two Euler equations. Assuming that the r-mode 
velocity field takes the form 

< = Ke; + vgel + v^el,] e^("*+™'^) , (102) 

and that u;,\g is similar, although of higher slow-rotation order than w,\ [32], we get 

iujpr^ sin^ 9v^ + imSp + 2pVlr^ sin 6* cos 9ve = , (103) 

and 

2iawr^ sin^ Bw^ + ims5T = . (104) 

The second relation shows iminediately that 5T is of higher slow-rotation order for the classic r-mode (since Wip is of 
higher order and w ^ fi for inertial inodes). This means that we can justifiably neglect the temperature variation in 
the pressure term in the first equation. That is, we have 

iojpr^ sin^ 9v^ -f- imSpo + 2/9J7r^ sin 9 cos 9vg = . (105) 

This is exactly the result used in J40]. Moreover, we now see that the first term in (|100p should also be of higher 
order, which means that we are left with 

V,<, « V,< . (106) 

The above relations allow us to estimate the terms needed to evaluate the bulk viscosity integrals once we have the 
leading order r-mode solution. Since the required contribution to the r-mode velocity field is purely toroidal, we have 

^« = -^^d^Yr , v^ = J^deVr . (107) 

r"' sm 6* r^ sm 9 

Moreover, the mode-solution is such that the only contribution comes from / = m in which case [/„ = (r/i?)™"'"^ 
(noting that the normalisation is irrelevant at the linear perturbation level). Finally, to the accuracy needed, the 
mode frequency is (in the rotating frame) 

- = ^ . (108) 

TTl + 1 

Using these results, and expanding the (scalar) pressure perturbation in spherical harmonics, i.e. using Sp^ = 
J2i SpiY{^e'^'^*+""f\ we find that ([1051) leads to 

^^.^^^C/_ (109) 

p V 2to 4- 3 m + 1 

and thus we have 

v,< « v,< « - ,^ ^, , ^1,. . ^^yr+i ■ (110) 

\/2m + 3 (rn + lj'^CQ 

This completes the approximation scheme, which provides all the expressions we need to make a rough estimate of 
the bulk viscosity damping timescale. 
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VII. MUTUAL FRICTION 



111 order to estimate the mutual friction damping, we need to determine the detailed counter-moving degree of 
freedom associated with the r-modes. In principle, this calculation proceeds as in [S^] and involves the rotational 
corrections to the shape of the star. However, in order to be consistent we will instead build on the approximations 
introduced in the previous section. First of all, it is easy to show that the countermoving degree of freedom will be 
poloidal (in contrast to the toroidal leading order r-mode solution). This means that we have 



1 



-WiY^ 



V, 



deYP 



Vi 



-.d^Yl' 



(111) 



If we also expand the temperature perturbation in such a way that 5T = ^j TiYf^ . then it follows from (|104l) that 



Vi 



2aijj 



Ti 



Similarly, we get from the radial component of (|30 



(112) 



W, 



2au} 



-drT, 



Combining these results with (|108p and (jllOp . we find that 



sT,n+l — 



SmaO,^ 



V277i + 3(m + l)3(2m + 5) 4 ^^^ 



R^ / r N^+i 



1 



r y 

'rj 



(113) 



(114) 



From this result we easily obtain the required counter-moving components of the r-modc. 

Wc arrive at (|114p by assuming that w^ = at the surface (in accordance with the discussion in [32|)- This condition 
is unlikely to be correct in the present problem, and a detailed analysis of the appropriate condition to impose on 
the phonons needs to be carried out in the future (c.f. the discussion in [2l|). At this point we simply note that our 
assumption that the whole star is composed of CFL matter is artificial anyway. Moreover, the boundary condition 
does not affect the parameter scaling of the result. This is the dominant factor in determining the r-mode damping 
timescale. 



VIII. THE R-MODE INSTABILITY WINDOW 



To estimate the growth/damping timescales for the r-mode instability we need to evaluate various energy integrals, 
see [2J, llol . First of all, the mode energy is given by 



1 r^ 

E = -m{m + l)i?-(2™+2) / pr2™+2rfr ^ 

2 Jo 

while the growth time-scale due to (current multipole) gravitational-wave emission is obtained from 

1 32^G172'»+2 (to-1)2'" /m + 2\""+' /-^ ,„,,,, 



c2m+3 [(2m-hl)!!]2 \m + l 



(115) 



(116) 



Turning to the shear viscosity, we know that |iy„s| <C l^n] for the superfiuid r-mode. This means that the shear 
viscosity integral simplifies to the usual, single fluid result. Hence, the leading-order energy loss due to shear viscosity 
follows from the expression 



Es- 



-m(m + 1) 




dr 



U„ 



+ {m-l)im + 2)\U,n\' 



dr 



'm{m^ - l)(2?7i + l)i?-(2'™+2) / jjr^"^dr , (117) 

Jo 



leading to 



1 



-(7n-l)(2m+l) / r]r'^"'dr 
Jo 



pr'^rn+2^^ 



(118) 
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Next, we need to estimate the bulk viscosity damping. From. (|44p we know that the general bulk viscosity expression 
is complicated, especially since we cannot rule out the possibility that the three different bulk viscosity contributions 
arc of similar magnitude. However, the approximations that we have made simplify the problem considerably. Ne- 
glecting the terms that arise due to the change in shape of the rotating star [the last line in (jH])], and making use of 
(fT06)) we find that 

iibv- J^ C,cs\yrv^\ av ^ (m + l)4(2m + 3)^ X 4 ^' ^^' 



which leads to 



_L _ 16m»4 [^ <5£,2™+4^^ 



Tbv (m + l)5(2m + 3) 7o % 



R 




(120) 



where 



Ccff « p'Ca • (121) 



This result follows after some surprising cancellations, associated with (|106p . It is notable that only one of the three 
bulk viscosity coefficients plays a role in our simplified case. Moreover, it is not the one {(2) that remains in the 
Navier-Stokes limit. Of course, one would not expect this drastic simplification for a more realistic stellar model. 
This fact provides strong motivation for a more detailed, numerical, r-mode analysis. We also learn that we can not, 
in general, ignore the additional bulk viscosities in the multi-fluid problem. 

Finally, we have the dissipation integral for the mutual friction. From the results discussed in |3ll.l32|. we find that 



E^f « -2 / pBn [sr - n"'n,j wi,,w';^dv . (122) 

Working out the angular integrals, this reduces to 

Emf « -2 1 pBil { [1 - g,^„+2 - Ql^,] Wl+, + [(m + l)(m + 2) - {m + 1)'QI+2 - (™ + 2)'QI+,] F^+J dr , 

(123) 
where we have used 

2 _ {l + m){l-m) 
^'~ {2l-l){2l + l)- ^^^^' 

The integral in (|123p is easily evaluated using the eigenfunctions from the previous section. 

In order to evaluate the relevant integrals, and obtain estimated timescales with explicit parameter scaling, we now 
assume that the density is uniform. Taking all parameters constant, we easily find that 

_/l.4Af0\ /10km\V P V / <. ^ 

for the I ^ m = 2 r-mode (the sign indicates that the mode is unstable). We also get, for the shear viscosity. 



. 1.1 X lO-^'' ^ — — ^^^^^^ s-^ , (126) 



and for the bulk viscosity we find 



5 /, . , r \ /-, \ 4 



1 _oo., in-3U Ccff \( R \ / l-47\/0\ / 1 ms ^ _i 



.. - '-' x '°-^' [TTJ^.) iTSTiS^J ItTJ [^ I -' ■ C^" 



Finally, in the case of the mutual friction we find 



^ « 0.176 f—^^ (-^] s-V (128) 

T„,f VlOknW Vims' ^ ' 
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FIG. 1: The r-mode instability window for a dense core comprising a CFL superconducting quark condensate and (finite 
temperature) phonons. The figure shows the critical rotational period vs core temperature for the instabiUty. Gravitational 
waves drive the I — m — 2 r-mode unstable for systems located below the different curves. The effect of shear- and bulk 
viscosity associated with phonon processes are shown as blue (dashed) and red (solid/dashed) curves, respectively. We show 
the result for both the canonical value /3 = 1 (solid red line) and the rather extreme value /3 = 10~^ (dashed red line) in (|86[) . 
leading to relaxation times of 10^ s and 1 s (as indicated in the figure), respectively. We also show the result for the case when 
a is increased by a factor of 10* (thin solid red line). Our result for the mutual friction damping (this dash-dot black line) is 
compared to the result from [la | (thin dash-dot black line) , demonstrating the dramatic effect of the different scaling with the 
rotation. For comparison, we also provide the instability curve for the combined shear- and bulk viscosity, obtained from (|87p 
and (|88[l . in the case of unpaired quark matter (black line at bottom of figure). Finally, we show the observed spin period and 
inferred core temperature for a number of astrophysical systems (filled squares), and we also indicate the region of parameter 
space where accreting neutron stars in low-mass X-ray binaries (LMXB) are thought to be located. The region below 10^'^ K, 
where the phonon-fiuid model may not be appropriate, is shown with grey background to emphasize the need for improved 
modelling in the astrophysically relevant part of parameter space. 



This result differs significantly from the estimate in [13|, in particular, in terms of the scaling with the rotation rate. 
However, our result is obtained from the actual counter-moving part of the r-mode solution and the final scaling 
accords with the standard result (sec, for example, |32|). The main difference between our result and that in [13| can 
be explained by noting that, if we take the r-mode energy to be proportional to w'^^^ rather than w^ then we arrive at a 
mutual friction damping time scale very similar to that obtained in [13| . The large discrepancy between the estimates 
thus follows from the fact that the velocites f ^ and luj^g associated with the r-modc solution arc of different orders of 
n in the slow-rotation approximation. 

Combining these estimates with ([75)) . ([77| . (|M)) and ([M|) . i.e. balancing the four timescales (|125p ~ (|128p . we arrive 
at the r-mode instability window in Figure [T] The figure shows the critical rotational period for the instability as 
function of the core temperature for a canonical star with mass 1.4 Mq and radius 10 km. Because of the less efficient 
damping mechanisms, we show the period rather than the rotation rate. This makes sense because the long timescales 
involved beome much clearer than when the critical rotation rate is expressed as a small fraction of the rotational 
break-up limit. As a result, gravitational waves drive the / = m = 2 r-mode unstable for systems located below the 
different curves in the figure. Our estimates show that the bulk viscosity is much weaker than the shear viscosity at 
all temperatures of interest. In order to show how the bulk viscosity result depends on the relaxation time, c.f. ()86p . 
we show the result for both the canonical value {3 — 1 and the rather extreme value /? = 10^^, leading to relaxation 
times of 10^ s and 1 s, respectively. The fact that the associated instability curves do not differ much shows that 
the results are not very sensitive to this parameter. This is due to the strong scaling with temperature. In order to 
explore the importance of the overall strength of the bulk viscosity, we also show the result for the case when a is 
increased by a factor of 10^ compared to our canonical value a = 1. This result shows that the bulk viscosity has to 
be vastly different from our assumed model in order to be relevant. The data illustrated in figure [1] also shows that 
the instability curve obtained from balancing the gravitational-wave driving and our result for the phonon mutual 
friction damping is vastly different from the result in [l3|. The different dependence on rotation is obvious, and it is 
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clear that the phonon mutual friction would be irrelevant for all astrophysical stars with CFL cores. For comparison, 
we show the instability curve for the combined shear- and bulk viscosity, obtained using ([57)) and (|55| , in the case of 
unpaired quark matter. It is interesting to note that the phonon shear viscosity in the CFL case leads to the strongest 
r-mode damping at temperatures above 10"'^*' K. Of course, astrophysical neutron stars will cool to temperatures much 
lower than this soon after birth so this region may not be that relevant. The results in the figure show, quite clearly, 
that we need to improve our understanding of the low-temperature regime. The instability curve below the 10^*^ K is 
obtained using our phenomenological model (j77p . Future work needs to model this regime in detail. 

For comparison, we also show the observed spin and inferred core temperature for a number of astrophysical systems 
in Figure [1] (the data is taken from [43] ) . This comparison is somewhat inconsistent because we have assumed that 
the core is shielded by a normal crust with a thickness of, at least, a few hundred meters (so that the standard 
heat-blanket argument |43| applies. The presence of a normal matter outer region should provide addition r-mode 
dissipation channels that we have not considered. In particular, one would expect a viscous boundary layer at the 
crust-core interface to be a more efficient damping agent (see [4J, [4^ for the most recent discussion) than the phonon 
mechanisms that we have considered. We also indicate the region of parameter space where accreting neutron stars in 
low-mass X-ray binaries would be located (given the usual arguments [4^, |43])- These systems would be at variance 
with the pure CFL plus phonon model in the sense that they would be located deep into the unstable region. If there 
were an active r-mode instability in these systems, one would expect it to have observational effects on, for example, 
the spin evolution. Observations do not provide any evidence of this. 

IX. CONCLUDING REMARKS 

We have presented the first true multi-fluid analysis of a dense neutron star core with a deconfined, colour-flavour- 
locked superconducting, quark core. By focussing on a cool system, and accounting only for the condensate and 
(finite temperature) phonons, we made progress by taking over much of the formalism from the analogous problem 
for superfluid ^He, the archetypal two-fluid laboratory system. The additional fluid degree of freedom, in the present 
case represented by the phonon gas, leads to the system not being well represented by the Navier-Stokes equations. 
In particular, a complete model requires a number of additional viscosity coefficients. 

Without an actual calculation it is not easy to establish whether the multi-fluid aspects are relevant or not. For 
example, in the case of the gravitational-wave driven instability of the f-modes it is known that the superfiuid degree 
of freedom is very important, since the vortex mutual friction may completely suppress that instability below a critical 
temperature (see [3l|, |43|). It is known that, because of the different nature of the associated velocity field, the mutual 
friction does not affect the r-mode instability in the same drastic fashion [S^, 143 • These examples provide clear 
evidence that different problems need to be considered on a case-by-case basis. 

We have provided a detailed dissipative formulation for a system comprising a quark (CFL) condensate and phonons. 
The model builds on recent improvements in our understanding of the analogous problems of superfluid ^He [9| 



and causal heat conductivity [21| . A key ingredient is the massless entropy component that represents the phonon 



excitations. We have discussed how the superfluid constraint of irrotationality reduces the number of required viscosity 
coefficients to four (one shear and three bulk) , and provided a translation between results in the literature [ll|, [i2| 
and the coefficients in our formalism. We also emphasised that many more dissipative channels may come into play 
in a rotating system where superfluid vortices are present |9l. Il9|. In particular, we showed how the vortex mediated 
mutual friction is accounted for in the model, and translated the available results for the associated coefficients [l3j . 

In order to be able to make relevant estimates for the r-mode instability, we developed a simple two-component 
equation of state based on the MIT bag model at zero temperatures with an additional phonon gas representing the 
thermal component. This example highlights the additional information that is needed in a multi-fluid analysis, in 
particular, regarding the entrainment coupling. Future work needs to provide a consistent equation of state, including 
all the key aspects. The model equation of state completed the formulation of the problem and we could, in principle, 
have carried out a numerical study of the r-modes. We opted not to do this, instead introducing a sequence of 
simplifying assumptions, because we felt that it would be useful to start by working out some less precise estimates 
for the relevant dissipation timescales. A numerical analysis of the problem should, of course, be encouraged. The 
problem has a number of interesting aspects, and may shed light on how one should deal with finite temperature 
superfluid neutron stars in general. 

This work was motivated by a desire to understand the different phases of CFL matter from a hydrodynamics 
point of view. It is, obviously, an interesting problem and the notion that observations of gravitational waves from 
relativistic stars may help shed light on the extreme QCD phase dia gram is exciting. The r-mode instability has been 
discussed in this context for some time, see for example |lll - [l3L [isl . [la . IstI . [SOI ] , but there has not been any previous 
discussion of the multi-fluid aspects of the problem. We hope that this work will stimulate a more detailed discussion 
between experts in the relevant areas. 
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The fact that the various phonon processes that we have accounted for can be shown to have httle effect on the 
r-mode instabihty should not discourage future efforts. In fact, the result could probably have been anticipated. The 
real challenge will be to account for additional degrees of freedom, especially associated with the kaons (either thermal 
or in a condensate) [15lll6l|. and the modelling of "hybrid" stars with quark cores and various phase transitions. These 
problems have additional features that, while possibly understood in principle, have never been considered in practice. 
This makes the modelling more complex, but also intriguing since the richer dynamics may lead to surprises. 

Appendix: Translation to the orthodox framework 

The relationship between the flux conservative formulation and the orthodox framework for superfluid Helium has 
been examined in detail elsewhere [§, [la| • Still, it is useful to summarise the main results that are needed to relate 
the different dissipation coefficients. First of all, we identify the "normal" fluid in the standard description with the 
gas of excitations, which is directly associated with the entropy of the system. This leads to 

Vn=<- (129) 

Secondly, the so-called "superfluid" velocity is given by 

vh^^^il-sK + evl, (130) 

where e = — is the entropy entrainment parameter. This means that the total mass flux takes the form 

p< = p^v'pf + psvs , (131) 

from which we learn that 

PS = -; and p^ = , (132) 

1 — £ 1 — £ 

where ps and pN sue the superfluid and normal fluid densities. 
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